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ABSTRACT 

We present Very Long Baseline Array observations of the kinematics of the 
water masers associated with OH 12.8—0.9, the fourth member of the so-called 
"water fountain" class of sources. We find that the masers occupy two distinct 
regions at the ends of a bipolar jet-like structure oriented north-south, with the 
blue-shifted masers located to the north and the red-shifted masers to the south. 
The masers are distributed along arc-like structures 12-20 mas across oriented 
perpendicular to the separation axis with an angular separation of ~110mas on 
the sky. Our multi-epoch observations, show the two maser arcs to be expanding 
away from each other along the axis of separation. The relative proper motions 
of the two maser regions is 2.7 mas yr^^ (~105kms~^ at the assumed distance 
of 8 kpc). The measured radial velocity difference between the northern, blue- 
shifted masers and the southern, red-shifted masers is 48.4 km s~^. The radial 
velocity, when combined with the proper motion, yields a three-dimensional ex- 
pansion velocity of 58kms~^ and an inclination angle of 24° for the jet. By 
combining our radial velocities with historical values, we estimate the three di- 
mensional acceleration of the masers to be ~0.63kms'^ yr~^ and a dynamical 
age for the coUimated outflow of ~90 yr. 



Subject headings: circumstellar matter — masers — stars: AGB and post AGB 
— stars: mass-loss — stars: individual (OH 12.8—0.9) 
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INTRODUCTION 



In the study of stars and stellar evolution, certain stages in the stellar life cycle remain 
a mystery simply because of their transitory nature. The evolutionary stage between the 
end of the asymptotic giant branch (AGB) and the planetary nebula (PN) phase is one such 
example. As stars evolve up the AGB, they lose mass at an increasing rate and spherically 
symmetric circumstellar envelopes (CSEs) are formed. Despite the spherical symmetry of 
the mass-loss process, a large fraction (~75%) of p lanetary nebulae (PNs ) show aspherical 
(i.e. elliptical, bipolar, quadrupolar) morphologies (IManchado et al.l 120001 ). The shaping of 
asymmetrical PNs must therefore occur over the relatively short time period of 10^-10'^ y r 



that the star spends in a post- AGB or proto-planetary nebul a (PPN) phase (IKwok 



Based on comparative studies of compact and extended PNs, Aaquist fc Kwokl ( 1996 ) find 



993|). 



that the lack of difference in the morphologies between the two groups s uggests that PNs 
morph ologies are primarily inherited from the AGB evolutionary stage. ISahai fc Trauger 
( 119981 ) surveyed a number of young PNs with the Hubble Space Telescope Wide Field Plan- 
etary Camera 2 and found that the majority of them were characterized by multi-polar 
bubbles distributed in a point-symmetric fashion about the central star. They also found 
collimated radial structures and bright equatorial struc tures indicating the pres ence of jets 
and disks/torii in some objects. These observations led ISahai fc Traugerl (119981 ) to propose 
a mechanism in which high-speed collimated outflows operate during the late-AGB or early 
PPN stages of stellar evolution. During this phase, the jet-like outflows carve out an imprint 
in the spherical AGB wind, and it is this imprint that provides the morphological signature 
necessary for the development of aspherical PNs. 

There is a small number of evolved objects that have been shown to exhibit the types of 
collimated jet-like outflows described above. These objects, dubbed "water fountain" neb- 
ulae, are thought to be stars entering the post-AGB or PPN evolutionary stage. They 
exhibit both H2O and OH maser emission, however, the relative characteristics of the 
two maser species differ from those of the typical AGB star. These differences are ap- 
parent in both the spectral profiles of the two types of masers and in the unique spatial 
morphologies of the masers as measu red by radio interferometry. Prior to our study of 
OH 12.8—0.9 (IBoboltz &: Maryell 12005), there were thr ee confirmed water-f ountain sources: 
IRAS 16342 -3814 JSahai et al.lll999l : lMorris et allboosh . IRAS 19134+2131 dlmai et al.l[20ol 
and W43A Jimai et al.ll2002h . 



Spectrally, these water-fountain sources exhibit double-peaked OH and H2O maser pro- 
files with the peaks symmetrically distributed about the radial velocity of the star. However, 
unlike the typical evolved star, the H2O maser peaks have a greater spread in velocity than 
the corresponding OH masers. For IRAS 16342-3814, IRAS 19134+2131, and W43A these 
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velocity ranges for the H2O masers are 259kms~^, 132kms~^ and 180kms~^ respectively 
(ILikkel et al.lll992l ). In addition, the expansion velocities (>60kms~^) implied by the maser 
spectra are much greater than the expansion velocity for the circumstellar wind of a typical 
AGB star (5-30 km s"!). 

Spatially, the H2O masers associated with these sources are known to trace bipolar 
jet-like outflows, hence the "water fountain" name. In the first Very Long B aseline Inter- 
ferometry (VLBI) study of the H2O masers toward W43A, llmai et al.l (120021 ) showed that 
the water masers are formed in a collimated precessing jet with an e stimated 3-D out 
fiow velocity of 145 k m s~^. Interferometric st udies of IRAS 163 42—3814 (IMorris et al.l 12003 



Claussen et al.l l2004f ) and IRAS 19134+2131 fllmai et al.l l2004f ) also found the H2O masers 
to exhibit bipolar distributions. Three-dimensional outfiow velocities for IRAS 16342—3814 
and IRAS 19134+2131 are 185kms~^ and 130kms~^, respectively. The bipolar jets traced 
by the H2O masers, presumably along the polar axis of the star, may represent the on- 
set of the axisymmetric morphologies that typify PNs. The dynamical ages of the jets for 
IRAS 19134+2131 and W43A are estimated to be ~50 yr and ~40 yr, r espectively (llmai et al. 



2OO2I . I2OOJ). That for IRAS 16342-3814 is estimated to be ~150yr (IClaussen et al.ll2004J ). 



The enigmatic source OH 12.8—0.9 was first classified by iBaud et al.l (119791 ) as a Type 
II OH/IR star based on its characteristic do uble-peaked 1612-M Hz OH maser profile. The 
distance to the source is unknown although iBaud et al.l (119851 ) included it in a sample of 
OH/IR stars associated with the Galactic center. OH 12.8— 0.9 has been linked to the in - 
frared source IRAS 18139— 1816, which is ~26arcsec away (Ite Lintel Hekkert et al.lll989l ). 
The SIMBAD astronomical database, however, still lists two separate entries for the source, 
one for each identifier. Classification of IRAS 18139— 1816 based on IRAS Low Resolution 



Spectrometer (8-23 fim) data was performed by iKwok et al. 



(1997) who placed the source 



Kwok et al.l ( 



19971 ) also noted 



in category /, a group with noisy or incomplete spectra, 
that IRAS 18139— 1816 is a "25 /xm peaker", or a source whose IRAS fiux density at 25 /im 
is greater than its fiux density at both 12 and 60 //m. Other 25 /xm peakers in clude young 



plane tary nebulae and carbon stars with circumstellar silicate emission features (iKwok et al. 



I997I ). OH 12.8—0.9 has been also observed with the PHT-S (2.5-11.6 n m) spectroph o tome- 
ter onboard the Infrared Space Observatory (ISO). From the spectra, iHodge et al.l (120041 ) 
classify the source in the 4/5- S A category. Sources in the 4/5 group peak long- ward of the 
11.6 /im limit of the PHT-S. Objects in the 4/5-SA sub-group additionally exhibit a deep 
10 /im absorption fea ture and absorption features from H2O at 3 and 6 //m, CO2 at ~4.3 //m, 
and CO at ~4.6 //m jHodee et al.ll2004h . 



H2O maser emission from OH 12.8—0.9 was first detected by lEngels et al.l (119861 ) who 
noted the fact that the H2O emission peaks are outside the range of the OH maser emission. 
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Gomez et al.l (119941 ) used the Very Large Array (VLA) to show that the "anomalous" H2O 
and OH mas er emission was spat ially coincident to within 1", and therefore belonged to the 



same source. 



Gomez et al.l (Il994l ) also found double-peaked profiles for both the OH and the 
H2O with the peaks of the OH separated by ~23kms~^ and H2O peaks separated by nearly 
twice this amount, ~42kms~^. Althoug h the v elocity range for the H2O is not as wide as 



the other water-fountain sources, lEngelj (120021 ) found that the shape and variations in the 
spectra were consistent with the water-fountain class and compatible with an axisymmetric 
wind. 



In lBoboltz Sz Marvell (120051 ) we reported on our initial Very Long Baseline Array (VLB A) 
observations of the H2O masers associated with OH 12.8—0.9. We found that the H2O masers 
trace the bipolar morphology typical of the water-fountain class. This bipolar structure com- 
bined with the spectral characteristics of the H2O and OH maser emission led us to propose 
that OH 12.8—0.9 is a fourth member of the rare water-fountain class of objects. In this 
article, we present additional multi-epoch observations made with the VLBA with the goal 
of measuring the kinematics of the H2O masers that trace the jet of OH 12.8—0.9. 



2. OBSERVATIONS 

We observed the 22.2 GHz H2O maser emission from OH 12.8-0.9 (a = 18''16™49^23, 6 = 
-18°15'01".8, J2000) using the 10 stations of the VLBA. The VLBA is operated by the Na- 
tional Radio Astronomy Observatory (NRAO)Q. VLBA spectral-line observations occurred 
on 2004 June 21, 2005 July 21, 2005 Nov. 2, and 2006 Feb. 24. Each observing run was 
approximately 5 hrs in length. A reference frequency of 22.23508 GHz was used for the 
H2O maser transition. Data were recorded in dual circular polarization using two 8-MHz 
(112.6kms~^) bands centered on the local standard of rest (LSR) velocity of — 58.0kms~^. 

For an unknown reason, the observations made on 2005 July 21 did not detect the H2O 
masers from OH 12.8—0.9. Following this non-detection, we requested a short (5 minute) 
observation of OH 12.8—0.9 with the Green Bank Telescope (GBT) in order to decide whether 
the masers would still be detectable in the pending two epochs of VLBA observations. On 
2005 Sept. 13 the GBT observed the H2O maser spectrum and found masers at levels of 2- 
3 Jy. We therefore continued with remaining two scheduled VLBA epochs and subsequently 
managed to detect and image the masers. The remainder of this article relates only to the 
three epochs in which the maser observations were successful, hereafter denoted: epoch 1 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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(2004 June 21), epoch 2 (2005 Nov. 2), and epoch 3 (2006 Feb. 24). 

The data were correlated at the VLBA correlator operated by NRAO in Socorro, New 
Mexico. Auto and cross- correlation spectra consisting of 512 channels with channel spacings 
of 15.63 kHz (~0.22kms^^) were produced by the correlator. Calibration of each of the three 
epochs was performed in accordance with standard VLBA spectral-line procedures using 
the Astronomical Image P rocessing System (AIPS) maintained by NRAO. The calibration 



of epoch 1 is described in iBoboltz fc Marvell (120051 ) and the two subsequent epochs were 



calibrated in a similar manner. For each epoch, residual delays due to the instrumentation 
were corrected by performing a fringe fit on the continuum calibrator (J1751+0939) scans. 
Residual group delays for each antenna were determined and applied to the target data. 

The bandpass response was determined from scans on J1751-I-0939 and was used to cor- 
rect the OH 12.8—0.9 data. The time-dependent gains of all antennas relative to a reference 
antenna were determined by fitting a total-power spectrum (from the reference antenna with 
the target source at a high elevation) to the total power spectrum of each antenna. The abso- 
lute flux density scale was established by scaling these gains by the system temperature and 
gain of the reference antenna. Errors in the gain and pointing of the reference antenna and 
atmospheric opacity variations contribute to the error in the absolute amplitude calibration, 
which is accurate to about 15-20%. 

Residual fringe rates were obtained by fringe-fitting a strong reference feature in the 
spectrum of OH 12.8—0.9. For each epoch we used the same strong feature at channel 
velocity Vlsr = —81.6 kms~^. The resulting fringe-rate solutions were applied to all channels 
in the spectrum. An iterative self-calibration and imaging procedure was then performed to 
map the emission in the reference channel for each epoch. The resulting residual phase and 
amplitude corrections from the reference channel were then applied to all channels in the 
8-MHz band. All of the above calibrations were then applied to the data prior to imaging. 



Imaging of the epoch 1 data is discussed in iBoboltz &: Marvell (120051 ) . and a similar 
methodology was employed for the two subsequent epochs. For epochs 2 and 3, full resolution 
images 2048 x 2048 pixels (~160x 160 mas) were generated using a synthesized beam sizes 
of 1.08 X 0.39 mas for epoch 2 and 0.88 x 0.33 mas for epoch 3. Images were produced 
for all spectral channels from — 88.3kms~^ to — 25.1kms~^ forming an image cube of 301 
image planes. Since peak maser flux densities were on the order of 1 Jy, channel images 
are noise limited and do not suffer from dynamic range limitations as is often the case in 
maser observations. Typical la off-source noise estimates for individual image planes were 8- 
10 mJy beam~^ for epoch 2 and 10-12 mJy beam~^ for epoch 3. The analysis and extraction 
of relevant information from the image cubes is described below. 



- 6 - 



3. RESULTS 

In order to identify and extract maser component parameters, two-dimensional Gaussian 
functions were fit to tlie emission in tlie image plane of each spectral (velocity) channel using 
the AIPS task SAD. Image quality was assessed using the off-source rms noise in the image. 
A cutoff flux density was set at five times the rms noise in the image plane containing the 
emission. Emission features with flux densities greater than the cutoff were fit with Gaussians 
to determine maser component parameters. The errors in right ascension and declination 
of the Gaussian fits to the identi fied componen ts were computed by the AIPS task SAD 



following the method outlined by I Condon! (119971 ). These position errors ranged from 5 /las 



for features with high signal-to-noise, to 110 /ias for features with lower signal-to- noise. 

Figures [1] and [2] show the spectral (upper sub-panels) and spatial (lower sub-panels) 
distributions of the H2O masers toward OH 12.8—0.9 from the analysis of our VLBA images 
for epoch 2 (Fig. [ Tl) and epoch 3 (Fig. [2]). A similar plot for epoch 1 was presented in 



Boboltz &: Marvell (120051 ). In the figures, panel (a) shows the entire range of H2O maser 
emission from the star. Panels (b) and (c) show enlarged views of the blue-shifted masers 
to the north and the red-shifted masers to the south, respectively. 

In the upper sub-panels in Figures[T]^a) and[2](a), the H2O maser spectra for epochs 2 and 
3 are plotted along with t he spectrum of the 1612-MHz OH masers from previous observations 
dSoboltz fc Marvell boosh as a reference. The velocities of the OH maser pe aks, have not 



chang ed significantly over time when compared to previous observations of iGomez et al. 



( I1994I ). We observed the velocities of the blue- and red-shifted OH peaks to be —68.0 and 



43. 7 km s ^, respectively with a vel ocity resoluti o n of 0.18 km s ^. By comparison, the 



velocities of the two peaks reported in lGomez et al.l (119941 ) are essentially the same at —67.4 



and — 44.0kms~^, to within the velocity resolution of their VLA observations (l.lkms^^). 
As in epoch 1, the H2O masers form a double-peaked profile with a velocity extent greater 
than that of the double-peaked OH maser profile. On the blue-shifted side of the emission, 
the peak flux density of epoch 2 remained the same as for epoch 1 at 1.3 Jy beam~^ and 
decreased slightly in epoch 3 to 0.97 Jy beam~^. For the red-shifted emission, the peak flux 
density steadily decreased from 1.0 Jy beam~^ in epoch 1 to 0.6 Jy beam~^ in epoch 2 to 

0. 2 Jy beam^^ in epoch 3. The velocity distributions for the masers in epochs 2 and 3 are 
very similar to those in epoch 1 and will be discussed in detail in § 13. 2[ 

The lower sub-panels of Figures [U^a) and [2]^a) show the entire spatial extent of the 
H2O masers for epochs 2 and 3, respectively. Enlarged views of the north and south regions 
are also shown in panels (b) and (c) of each figure. We find that in epochs 2 and 3 the 
overall spatial structure of the H2O masers has remained essentially unchanged from epoch 

1. There are two distinct regions oriented roughly north-south on the sky separated by 
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~110mas. The position angle of the axis of separation between the mean position centers 
of the northern and southern maser regions is ~2.6° east of north. The masers in the two 
regions are arranged in arc-hke structures oriented roughly east-west with the exception 
of the feature located near (15.0, —10.0) in Figures [T](b) and[2|^b). This feature, although 
apparent in the epoch 1 image c ube, did not meet all of the component selection criteria in 
our initial analysis of the data (IBoboltz fc Marvell l2005l ). The feature is discussed further 
in § 13.11 below. The blue- and red-shifted arcs of masers are approximately 20 and 12 mas 
across, respectively, corresponding opening angles for the jet of 13-20°. 



3.1. Maser Proper Motions 

In Figures [1] and [21 we plotted the maser emission identified in every spectral channel, 
thus an individual maser feature will appear in multiple spectral channels and will consist 
of multiple points in the spatial distributions shown in the figures. In order to track and 
estimate the motions of the maser features, it is necessary to determine average values in right 
ascension, declination, and velocity for each feature. These average values were computed 
using a flux-density-squared weighted average for emission identified in two or more channels 
and spatially coincident to within 1 mas. The flux assigned to the maser feature was simply 
the peak flux in the channels spanned by the feature. Characteristics of the masers for the 
three epochs are presented in Tables [TH5J 

Since we did not use the technique of phase-referencing for our three VLBA epochs, 
the absolute position of the phase center in each image cube is unknown. However, in the 
reduction and analysis of each epoch of observations, we used what we believe to be the 
same maser feature as a phase reference. This feature was the strongest feature in each of 
the three epochs. In addition, it maintained the same velocity and spatial location relative 
to the other three northern, blue-shifted masers common to all three epochs. The averaging 
of the maser identifications over several channels resulted in slight (<0.1 mas) position shifts 
for the reference feature in each of the three epochs. The coordinate frames for all three 
epochs were therefore realigned such that the reference feature (the third blue-shifted feature 
in Tables [1] and the second blue-shifted feature in Tables [2] and [3]) coincided with the origin. 

Relevant maser characteristics determined from the analysis are plotted in Figures [3] 
and m Figure [3] shows the blue-shifted H2O masers in the north of OH 12.8—0.9. In this 
figure, the alignment of the reference features for the three epochs at the origin is apparent. 
Also common to all three epochs are features near (—2.0, 0.5), (8.0, 0.2) and (15.0, —10.0). 



As me ntioned earlier, the feature near (15.0, —10.0) did not appear in iBoboltz fc Marvel 



(I2OO5I ) because it did not meet all of the selection criteria at the time. However, with the 
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additional knowledge from epochs 2 and 3, and some reprocessing of the epoch 1 image 
cube, we were able to extract the relevant parameters for this feature. From Figure [3] it is 
apparent that any position differences between the three epochs are small, thus indicating 
little motion of the northern features relative to the reference feature over the course of the 
three experiments. This, however, is not the case for the southern maser features. Figure H] 
shows the red-shifted H2O masers in the south of OH 12.8—0.9. Here there are five features 
that are common to all three epochs. Clearly there is a consistent proper motion of all 
of the features nearly due south relative to the blue-shifted reference feature. The average 
angular separations of the north and south masers are 107, 11 and 111 mas for epochs 1, 2 
and 3, respectively. As discussed in iBoboltz fc Marvell (120051 ) the distance to OH 12.8—0.9 
i s not well knowii . Best estimates place the object near the Galactic center at D ^ 8kpc 
( iBaud et al.lll985l ). Assuming this distance, then the linear separations between the blue- 
and red-shifted masers are approximately 860, 880 and 890 AU, respectively. 

In order to better characterize the net expansion of the masers we computed separations 



between pairwise combinations o 



characterize O H (IChapman et al. 



components. This technique has previousl y been used to 



1991 



Bloemhof et al.lll992l : lKemballlll992f ). H2O (IMarvel 



1996h . and SiO JBoboltz et al.lll997l " IChen et al.ll2006l ) maser motions and has no dependence 



on the alignment of the maps or a priori assumptions about the velocity field. The procedure 
involves computing the angular separation {A9) between two features at one epoch (epoch A) 
and the separation between the corresponding two features at a second epoch (epoch B). 
The difference between the two values of is the pairwise separation and can be written 

A6'b - A^A = - TjIb - In - TjIa , i = l,n; j = i + l,n (1) 

where = {xi,yi) and {xi,yi) are the relative offsets in right ascension and declination, 
respectively. The procedure is repeated for all possible pair combinations and the separations 
can be plotted as a histogram (e.g. Figure E]). The inclusion of all possible pair combinations 
often results in bimodal distribution with one of the peaks biased towards zero. This is 
because some of the separations involve pairs of closely spaced maser components on the same 
side of the distribution (i.e. the northern or southern regions for OH 12.8—0.9) that have 
little motion relative to one another. For the sake of clarity, and to determine representative 
values for the angular shifts due to the expansion, we have included only those pairs separated 
by more than 80 mas in the histograms shown in Figure [51 Figure |5] plots the pairwise 
separations over time intervals of (a) 613, (b) 500 and (c) 114 days, respectively. All three 
histograms have centroids which are biased toward positive values indicative of expansion. 
The mean (median) differences between the north and south masers are 4.64 (4.69) mas, 3.88 
(3.89) mas, and 0.77 (0.80) mas for 613, 500, and 114 days, respectively. The la standard 
deviations are 0.28, 0.21 and 0.13 mas, respectively. The equivalent mean angular velocities 
are 2.76 ± 0.16 mas yr~^, 2.83 ± 0.15 mas yr~^ and 2.47 ± 0.41 mas yr^^. To within the 
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errors, the three velocities are consistent, and it is impossible to determine whether there is 
any acceleration of the masers in the outflow. 

If 8 kpc is again assumed as the distance to OH 12.8—0.9, the linear proper motions of 
the masers can be computed for the angular velocities above. The resulting average linear 
separation velocities are 105 ± 6 km s~^, 107 ± 6 km s~^ and 94 ± 16 km s~^ for the 613, 
500, and 114 day intervals respectively. If we assume that the north and south masers are 
moving outward at equal speeds (i.e. at half the computed separation velocity), then the 
outflow velocity in the plane of the sky would be ~53 km s~^ for the 613 day interval. This 
velocity is slightly more than double the outflow velocity determined for OH 12.8—0.9 in 



Boboltz &: Marvell (120051 ) based solely on the radial velocities. In the following section we 
update the radial velocities of the masers and combine them with the proper motions in 
order to estimate the full 3-D kinematics of the masers in the jet. 



3.2. 3-D Maser Outflow 



In iBoboltz fc Marvell (120051 ) . we treated the radial velocity of OH 12.8— 0.9 by simply 
determining the velocity separation of the blue- and red-shifted maser features representing 
the peaks in the spectral distribution. This was an improvement over previous single-dish 
and VLA studies that contained limited spatial information and likely suffered from blending 
of features in the spectral domain. In light of the two new epochs of VLBA data, we felt that 
the more rigorous approach of computing spectral separations (Av) from multiple features 
was more appropriate. To do this we averaged the velocities of the blue- and red-shifted 
masers identified in Tables dHS] separately. We also computed the standard deviations for 
each set of components. From these averages we computed the spectral separation Av at 
each epoch. The error in A'u is simply the standard deviations added in quadrature. Using 
this method, we find mean values for Av of 48.7 ± 1.6 km s~^, 48.3 ± 1.3 km s~^ and 
48.2 ± 1.3 km s^^ for epochs 1, 2, and 3, respectively. These values are also listed in Table HI 



In iBoboltz &: Marvell (120051 ) we compared the spectral separation of the blue- and red- 
shifted peaks in our spectrum with peak separations fr o m prev i ous s ingle dish and VLA 



observations (i.e. lEngels et al.l Il986l : I Gomez et al.l 1 1994) : lEngeld |2002| ). Here we wish to 



compare our spectral s eparations com puted from the component averages with these previous 
results, however, only lEngelsl (120021 ) tabulates parameters for individual spectral features. 
Engeld (120021 ) categorizes components A-H of his Table 4 into the blue-shifted group and 
components M-Q of Table 5 into the red-shifted group, with the remaining components, I-L, 
listed as intermediate. From this information, we computed average velocities for the blue- 
and red-shifted spectral features and values for Af in a manner similar to that performed 
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for our own VLBA data. The resulting values of At; and the corresponding errors for the 
Engeld (|2002| ) are listed in Table |H The only drawback to the Engels data as compared to 
our VLBA data is the fact that there is no spatial information, thus the degree to which 
spectral blending is a factor is unknown. It should be noted, that the spectral resolution of 
the Engels measurements is slightly better than our VLBA resolution at 0.16 km s~^. 

Plotted in Figure [6] are the results of our computation of Ay frorn the bl ue- and red- 
shifted component averages. The values from the 11 epochs of lEngeld (120021 ) are plotted 
as filled triangles and those from our three epochs of VLBA data as filled circles. The 
data clearly shows an increase in At; as a function of time. To this data, we performed a 
linear least-squares fit weighted by squares of the errors. This fit is plotted as the line in 
Figure El with a slope of 0.53 ± 0.04 km s~^ yr~^ corresponding to the relative acceleration 
of the blue- and red-shifted masers in the radial direct i on. T his is slightly lower than the 
0.68 ± 0.06 km s"^ yr"^ reported in iBoboltz fc Marvell (|2005[ ). With the new procedure of 
averaging over multiple spectral features, we have attempted to remove the error involved 
in computing the velocity separation from single peaks on the blue- and red-shifted sides 
of the spectrum th at are l i kely u nrelated over over periods greater than the lifetime of the 
components < 3 yr lEngeld (120021 ). Engels also discusses ~1 km s~^ velocity shifts that were 
observed for some components in the spectra of OH12.8— 0.9 and the likelihood of these 
shifts being caused by intensity changes within blended features. Such shi fts could affect the 
velocity separations computed from single spectral peaks as was done in iBoboltz fc Marvel 
J2005h . 



If we assume an average over our three VLBA epochs, Av ~ 48 km s^^, as the speed 
the masers are moving away from each other along the line of sight to the observer, then 
the outflow velocity in the radial direction is ~24 km s^^. The 3-D outflow velocity of the 
masers may be estimated by combining the radial motion above with the motion in the 
plane of the sky determined in § 13.11 We find the 3-D outflow velocity of the masers to be 
~58 km s^^ with an inclination angle for the jet of ~24° with respect to the plane of the sky. 
The 3-D outflow acceleration, assuming this inclination angle and the radial acceleration 
above, is ~0.63 km s~^ yr~^- The outflow acceleration in the plane of the sky would be 
~0.58 km s~^ yr~^ under this assumption. Thus over the 613 days between epochs 1 and 
3, we could expect a change in velocity of the blue-shifted masers relative to the red-shifted 
masers of 1.93 km s~^ or 0.05 mas yr^^, which is undetectable in our present observations. 

From the above information and two additional assumptions, namely: that the dynam- 
ical center of the outflow is the midpoint along the axis of separation between the two maser 
regions and that the masers have zero initial velocity, then an upper limit to the dynamical 
age of the outflow can be computed. We find this dynamical age to be ~90 yr. This value 
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is roughly consistent with the age determined in iBoboltz &: Marvell (120051 ) , which was esti- 
mated with no knowledge of the maser motions in the plane of the sky or the inclination 
angle of the outflow. Assuming that the acceleration has and will remain constant at the 
above value of 0.63 km s~^ yr~^, then the total time to reach an outflow velocity comparable 
to the other water fountain sources {Av ~ 150 km s~^) is only about ~240 yr, roughly 2.5 
times the current dynamical age of OH 12.8—0.9. 

In order to compare the characteristics of OH 12.8—0.9 to the other water-fountain 
sources, we have summarized the relevant properties of all four sources in Table O From 
the table we see that the assumed distance to OH 12.8—0.9 is intermediate among t he four 



source s. The peak flux density i s shown as a r ange m easured over multiple epochs by lEngels 



f l2002[ ) for OH 12.8-0.9 and by iLikkel et all (119921 ) for the three other sources. The peak 



flux range for OH 12.8—0.9 is similar to that of the other distant source IRAS 19134+2131 
and weaker than those of W43A and IRAS 16342—3814. The angular extent of the outflow 
on the plane of the sky is also similar to that of IRAS 19134+2131. The linear extent of 
OH 12.8—0.9, however, is much less than that of any of the other water-fountain sources. 
An increase in the assumed distance to the source, would serve to bring this value in line 
with the other sources. For the coUimation of the OH 12.8—0.9 jet, we have used an average 
of the northern and southern arcs and find that the coUimation is slightly larger than for 
the other three water-fountain sources. If we use only the southern arc, or we disregard 
the one northern feature at (15.0, —10.0) then the coUimation is ~10°, which is consistent 
with the other water-fountain sources. The estimated inclination angle of the outflow is also 
comparable values estimated for the other three objects. 

Aside from the difference in linear extent, the primary difference between OH 12.8—0.9 
and the other water-fountain sources is in the kinematics. Both the radial (vrad) and tangen- 
tial (vtan) outflow velocitlcs are less than half the values for the other sources. It is possible 
that the tangential velocity could be closer to that of the other sources if the distance to 
the source turns out to be greater than 8 kpc. An increased distance would not increase the 
radial velocity estimate, but would instead result in a decreased estimate for the inclination 
angle of the jet. Another major difference between OH 12.8—0.9 and the three other water 
fountains is that OH 12.8—0.9 appears to be t he only objec t with a measurable acceleration 



Orad) of the masers along the line of sight. iLikkel et al.l (119921 ) observed the H2O maser 



spectra of the three other water-fountain sources over the course of a few years and found no 
evidence for systematic velocity drifts in spectral features that would suggest acceleration of 
the masers. The acceleration in combination with the smaller extent of the masers toward 
OH 12.8—0.9 might suggest that it is younger than the other three objects. This is not, 
however, confirmed by the computed dynamical age of the jet, which is intermediate among 
the four sources at 90 yr. 
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One characteristic that is not shown in the tables is the morphology of the masers rel- 
ative to the central star and the jet direction. From the images of OH 12.8—0.9 presented 
here, we find that the masers are arranged in point symmetric fashion along arcs perpen- 
dicular to the direction of the outflow, suggesting that the outflow is impacting a denser 
surrounding medium. This bow shock structure is also observe d for the H2O maser s toward 
IRAS 16342-3814 at least for the blue-shifted side of the jet Jciaussen et al.ll2004[ ). These 
two sources are in contrast to W43A where the masers are arranged along the direction of 
the outflow. The struct ure of the W43A masers is well represented by a precessing jet model 
(jimai et al.ll2002l . |2005| ). The remaining water-fountain source, IRAS 19134+2131, shows the 
point symmetric morphology of the three other sources, but the relationship between the 
masers and the o utflow is less clear. From the two epochs of VLBI observations reported in 
Imai et al.l (120041 ). the masers on the eastern red-shifted side appear to be roughly aligned 
with the line connecting them to the western blue-shifted masers. The motions of the masers, 
however, do not appear to be along this line, but rather perpendicular to this in a primarily 
north-south direction. 

Similar to the W43A outflow, precession is indicated for IRAS 16342—3814, but not as 
a result of the availa ble H2O maser da ta. Instead, it was the Keck adaptive optics images 
in the near-infrared (ISahai et al.ll2005l ) that showed a corkscrew structure inscribed on the 
walls of the observed lobes. The combined radio/infrared observations demonstrate that 
both bow shock morphologies and precession can take place in the same object. Since no 
such high-resolution optical/infrared imaging is available for OH 12.8—0.9 it is unknown 
whether precession occurs in the object, however, it is clearly not suggested by the structure 
or kinematics of the H2O masers. 



4. CONCLUSIONS 

Using the VLBA we have shown that, like the spatial morphology, the kinematics of the 
H2O masers toward OH 12.8—0.9 indicate that the object is a member of the rare "water 
fountain" class of sources. We find that the masers continue to be located in two distinct 
regions at the ends of a bipolar jet-like structure having a north-south orientation on the 
sky. The masers in the two regions are distributed along arc-like structures ~ 12-20 mas 
across oriented roughly perpendicular to the separation axis. The angular separation of the 
two regions is roughly 110 mas, with the blue-shifted masers located to the north and the 
red-shifted masers to the south. The arc-like arrangements of the masers are suggestive of 
bow shocks formed by a coUimated axisymmetric wind impinging on the ambient medium 
surrounding the star. 
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The two additional epochs of VLBA observations, beyond our first epoch in 2004, have 
allowed us to track the projected spatial motions of individual features over the course of 
~600 days. We find that the two maser regions maintain their arc-like appearance as they 
expand away from each other along the axis of separation. The observed proper motion 
of the southern masers relative to the northern masers is 2.76 mas yr~^ (~105kms^^ at 
the assumed distance of 8 kpc). The masers maintain coUimation during this expansion 
with little motion perpendicular to the north-south axis. Combining the proper motions 
with the measured radial velocities yields a three-dimensional (3-D) expansion velocity of 
58 km s~^ relative to the dynamical center of the distribution. The inchnation angle of the 
jet is estimated to be ~24°. A simple linear fit to the historical radial velocity data combined 
with our own measurements yields an estimate for the 3-D acceleration of ~0.63 km s~^ yr~^ 
and a corresponding dynamical age for the outflow of ~90 yr. Long-term VLBA monitoring 
of the H2O masers should enable the determination of the tangential acceleration of the 
masers and the true 3-D acceleration of the outflow associated with OH 12.8—0.9. 

This research has made use of the SIMBAD database, operated at CDS, Strasbourg, 
France. This research has made use of NASA's Astrophysics Data System Bibliographic 
Services. 
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Table 1. H2O maser characteristics derived from the epoch 1 (2004 June 21) VLB A data. 



^'LSR 

(kms~^) 


Si/ 

(Jy beam~^) 


(Jy beam~^) 


Relative R.A. 
(mas) 


ctr.a. 
(mas) 


Relative Dec. 
(mas) 


fDoc. 

(mas) 


J-Jl ollll. UCU 














—85.07 


0.044 


0.009 


5.869 


0.038 


2.520 


0.063 


—83.32 


0.052 


0.009 


14.806 


0.034 


— 10.133 


0.049 


-81.71 


1.297 


0.012 


0.000 


0.006 


0.000 


0.008 


-81.67 


0.123 


0.012 


-1.813 


0.031 


-0.322 


0.045 


-80.88 


0.261 


0.009 


7.647 


0.009 


0.376 


0.014 


Red-shifted 


masers 












-34.61 


0.247 


0.009 


1.116 


0.012 


-108.821 


0.019 


-33.92 


0.139 


0.012 


-5.773 


0.018 


-108.460 


0.027 


-33.72 


0.498 


0.012 


2.235 


0.018 


-108.856 


0.030 


-33.27 


1.041 


0.012 


-2.807 


0.005 


-108.597 


0.007 


-32.58 


0.499 


0.009 


5.297 


0.010 


-106.906 


0.017 



Table 2. H2O maser characteristics derived from the epoch 2 (2005 Nov. 2) VLBA data. 



VhSR 

(kms-^) 


(Jy beam~^) 


(Jy beam~^) 


Relative R.A. 
(mas) 


cr.a. 
(mas) 


Relative Dec. 
(mas) 


CDec. 

(mas) 


Blue-shifted 


masers 












-83.29 


0.153 


0.009 


15.458 


0.029 


-10.516 


0.045 


-81.57 


1.323 


0.010 


0.000 


0.017 


0.000 


0.028 


-81.41 


0.096 


0.010 


-2.207 


0.041 


-0.456 


0.063 


-80.73 


0.216 


0.010 


8.174 


0.032 


0.173 


0.051 


Red-shifted 


masers 












-34.47 


0.468 


0.010 


1.286 


0.020 


-112.898 


0.033 


-33.54 


0.293 


0.009 


2.518 


0.020 


-112.793 


0.034 


-33.31 


0.086 


0.009 


-6.432 


0.048 


-112.695 


0.071 


-33.10 


0.569 


0.009 


-2.658 


0.020 


-112.669 


0.030 


-32.68 


0.142 


0.009 


5.765 


0.034 


-110.732 


0.051 
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Tablc 3. H2O mascr characteristics derived from the epoch 3 (2006 Feb. 24) VLBA data. 



VLSR 

(kms-^) 


(Jy beam"^) 


(Jy bearn"^) 


Relative R.A. 
(mas) 


ctr.a. 
(mas) 


Relative Dec. 
(mas) 


CDec. 

(mas) 


Blue-shifted 


masers 












-83.29 


0.174 


0.013 


15.580 


0.029 


-10.741 


0.054 


-81.54 


0.970 


0.013 


0.000 


0.026 


0.000 


0.043 


-81.30 


0.079 


0.013 


-2.358 


0.063 


-0.384 


0.110 


-81.26 


0.065 


0.013 


0.347 


0.068 


1.034 


0.111 


-80.72 


0.094 


0.012 


8.208 


0.043 


0.067 


0.074 


Red-shifted 


masers 












-34.44 


0.228 


0.015 


1.315 


0.039 


-113.770 


0.072 


-33.46 


0.151 


0.014 


2.582 


0.039 


-113.587 


0.070 


-33.35 


0.061 


0.012 


-6.378 


0.062 


-113.447 


0.106 


-33.04 


0.164 


0.012 


-2.593 


0.041 


-113.582 


0.075 


-32.71 


0.123 


0.012 


5.717 


0.047 


-111.546 


0.085 
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Table 4. Velocity separation of blue- and red-shifted water masers. 



Obs. Date 


Av 


O'Av 


Spectral res. 




(yr) 


(kms^"'^) 


(kms""'^) 


(kms^""^) 


Ref. 


1987.452 


38.8 


1.5 


0.16 


1 


1988.516 


37.7 


1.6 


0.16 


1 


1990.132 


42.0 


2.1 


0.16 


1 


1990.249 


41.2 


2.3 


0.16 


1 


1991.332 


41.5 


2.1 


0.16 


1 


1992.049 


42.0 


2.0 


0.16 


1 


1992.975 


42.3 


2.3 


0.16 


1 


1994.184 


40.9 


2.1 


0.16 


1 


1995.186 


42.4 


2.0 


0.16 


1 


1996.675 


43.3 


2.3 


0.16 


1 


1999.088 


44.6 


2.0 


0.16 


1 


2004.473 


48.7 


1.6 


0.22 


2 


2005.838 


48.3 


1.3 


0.22 


2 


2006.151 


48.2 


1.3 


0.22 


2 



Note. — References: (1) 



Engels 



(|2002l ). (2) this work. 



- 19 - 



Table 5. Summary of the H2O maser characteristics for the water fountain sources. 



Characteristic 


OH 12.8-0.9 


W43A 


IRAS 19134+2131 


IRAS 16342-3814 


Distance (kpc) 


8 


2.6 


16 


2 


Speak (Jy) 


1.4-13.0 


11.5-56.6 


1.2-6.0 


16.0-67.7 


Angular Extent (mas) ^ . . 


110 


920 


135 


3000 


Linear Extent (AU) ^ 


880 


2400 


2180 


6000 


Outflow Collimation (deg.) 


15 


5 


10 


6 


Inclination Angle (deg.) 


24 


39 


25 


40 


Outflow Velocity: 










wtan (kms"^) 


53 


110 


120 




Wrad (knis^i) 


24 


90 


65 


130 


vtot (kms^^) 


58 


145 


130 




Outflow Acceleration: 










arad (kms^i yr^i) 


0.53 








Dynamical Age (yr) 


90 


50 


50 


150 ^ 


References 


1, 2, 3 


4, 5, 6 


4, 7 


4, 8, 9, 10 



'^Range over multiple epochs of single dish observations. 
''Measured in the plane of the sky. 
■^Relative to the plane of the sky. 
■^From OH maser observations. 



Not e. — References: (plEngelsl ( iooj). r2)lBoboltz fc Mar^j2005[) (3) this w ork. f4'llLikkel et"aL 
(Il992h. (5) llmai et al.l tooi) r6'l llmai et allHoO^ . (7) llmai et all (|2004h . (8) ISahai et al.l (Il999h . 
f9l lMorris et all |2003h . (10) Iciaussen et al.l 12004 1 
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Fig. 1. — The H2O maser emission toward OH 12.8-0.9 from epoch 2 (2005 Nov. 2). In panels 
(a), (b) and (c), the top sub-panels show the spectra formed by plotting maser component flux 
density versus LSR velocity, color-coded according to maser velocity. A dark solid line in each 
upper sub-panel represents the scalar-averaged cross-power spectrum on the Los Alamos-Pie Town 
VLBA baseline. A dotted line in the top sub-panel of (a) represents the scalar-averaged cross- 
power spectrum of the OH masers (with the flux density scaled by a factor of 0.5) from the Los 
Alamos-Pie Town baseline. The lower sub-panels in (a), (b) and (c) plot the spatial distribution of 
the H2O masers with point-color representing the corresponding velocity bin in the spectrum and 
point-size proportional to the logarithm of the maser flux density. Panel (a) shows all H2O maser 
components from our VLBA observations. A dashed line represents the axis of separation between 
the centers of the blue- and red-shifted masers at a position angle of 2.6° east of north. The "x" 
symbol represents the midpoint of the bipolar distribution. Panels (b) and (c) show expanded views 
of the blue- and red-shifted maser features respectively. Errors in the positions of the features are 
smaller than the data points for all panels. 
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Fig. 2. — Same as Figured], but for epoch 3 (2006 Feb. 24). 
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Fig. 3. — Northern blue-shifted masers associated with OH 12.8—0.9 for all three epochs 
(2004 June 21, 2005 Nov. 2, and 2006 Feb. 24). Masers are color-coded according to 
velocity and point-size is proportional to the logarithm of the maser flux density. There are 
four maser features common to all three epochs including the reference feature at (0,0). The 
masers show little motion relative to this stationary reference feature. 
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Fig. 4. — Southern red-shifted masers associated with OH 12.8—0.9 for all three epochs (2004 
June 21, 2005 Nov. 2, and 2006 Feb. 24) demonstrating the motions of the masers over time 
relative to the reference feature in Figure [31 Masers are again color-coded according to 
velocity and point-size is proportional to the logarithm of the maser flux density. All five 
features are common to all three epochs. 
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Fig. 5. — Histogram showing the pairwise spatial separations between combinations of north- 
ern and southern maser pairs for (a) epochs 3 and 1, (b) epochs 2 and 1, and c) epochs 3 
and 2. Listed in each panel is the elapsed time between epochs in days. 
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Fig. 6. — Plot showing the measured velocity separation {Av) of the blue- and red-shifted 
water maser features as a function of time from values listed in Table HI The line repre- 
sents a weighted linear least-squares fit to the data and indicates a constant acceleration of 
0.53 km s^^ JT^^. 



